Nonlinear electrodynamics and FRW cosmology 
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Maxwell electrodynamics, considered as a source of the classical Einstein field equations, leads 
to the singular isotropic Friedmann solutions. We show that this singular behavior does not occur 
for a class of nonlinear generalizations of the electromagnetic theory. A mathematical toy model is 
proposed for which the analytical nonsingular extension of FRW solutions is obtained. 
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I. INTRODUCTION 

The standard cosmological model, based on 
Friedmann-Robertson- Walker (FRW) geometry with 
Maxwell electrodynamics as its source, leads to a cosmo- 
logical singularity at a finite time in the past [Q. Such 
a mathematical singularity itself shows that, around the 
very beginning, the curvature and the energy density 
are arbitrarily large, thus being beyond the domain of 
applicability of the model. This difficulty raises also 
secondary problems, such as the horizon problem: the 
Universe seems to be too homogeneous over scales which 
approach its causally correlated region ||^. These sec- 
ondary problems are usually solved by introducing ge- 
ometric scalar fields (for a review on this approach see 
Ref. and references therein). 

There are many proposals of cosmological solutions 
without a primordial singularity. Such models are based 
on a variety of distinct mechanisms, such as cosmologi- 
cal constant Q , nonminimal couplings Q , nonlinear La- 
grangians involving quadratic terms in the curvature [ 3 1 , 
modifications of the geometric structure of spacetime [7|, 
and nonequilibrium thermodynamics j^, among others. 
Recently, an inhomogeneous and anisotropic nonsingular 
model for the early universe filled with a Born-lnfeld- 
type nonhnear electromagnetic field was presented 
Further investigations on regular cosmological solutions 
can be found in Ref. ||l0| . 

In this paper it is shown that homogeneous and 
isotropic nonsingular FRW solutions can be obtained by 
considering a toy model generalization of Maxwell elec- 
trodynamics, here presented as a local covariant and 
gauge-invariant Lagrangian which depends on the field 
invariants up to the second order, as a source of classical 
Einstein equations. This modification is expected to be 
relevant when the fields reach large values, as occurs in 



the primeval era of our Universe. Singularity theorems 
jm are circumvented by the appearance of a high (but 
nevertheless finite) negative pressure in the early phase 
of FRW geometry. In the Appendix we consider the in- 
fluence of other kinds of matter on the evolution of the 
universe. It is shown that standard matter, even in its 
ultrarelativistic state, is unable to modify the regularity 
of the obtained solution. 

Heaviside nonrationalized units are used. Latin indices 
run in the range (1,2,3) and Greek indices run in the 
range (0,1,2,3). The volumetric spatial average of an 
arbitrary quantity X for a given instant of time t is de- 
fined as 



a: lim - / X . 



(1) 



where V = J \/—g cPx'^, and Vo stands for the time de- 
pendent volume of the whole space. 



II. EINSTEIN-MAXWELL SINGULAR 
UNIVERSE 

Maxwell electrodynamics usually leads to singular uni- 
verse models. In a FRW framework, this is a direct conse- 
quence of the singularity theorems pl , and follows from 
the exam of the energy conservation law and Raychaud- 
huri equation p^. Let us set the line element 



ds'=c'dt' 



1 + er2/4c2 



(2) 

where e = —1, 0, -1-1 hold for the open, flat (or 
Euclidean) and closed cases, respectively. The 3- 
dimensional surface of homogeneity t = const is orthog- 
onal to a fundamental class of observers represented by 
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a four- velocity vector field = cS^^. For a perfect fluid 
with energy density p and pressure p, the two above- 
mentioned equations assume the form 



(3) 



III. NONSINGULAR FRW UNIVERSES 



The toy model generalization of Maxwell electro- 
magnetic Lagrangian will be considered up to sec- 
ond order terms in the field invariants F and G = 



-AE ■ H as 



-^4 k , ^ ^ 



(4) 



in which k is the Einstein gravitational constant and the 
overdot denotes Lie derivative respective to v, that is 
{\/c)d/dt. Equations ^ and (^) do admit a first inte- 
gral 



k A\ e 



(5) 



■F + aF^ 



(13) 



where a and (3 are arbitrary constants. Maxwell electro- 
dynamics can be formally obtained from Eq. ( p^ ) by set- 
ting a = = /3. Alternatively, it can also be dynamically 
obtained from the nonlinear theory in the limit of small 
fields. We will not consider generalizations of Eq. jl^ ) 
which include the term FG in order to preserve parity. 
The energy-momentum tensor for nonlinear electromag- 
netic theories reads 



Since the spatial sections of FRW geometry are 
isotropic, electromagnetic fields can generate such a uni- 
verse only if an averaging procedure is performed [ p^ . 
The standard way to do this is just to set for the electric 
Ei and magnetic Hi fields the following mean values: 



E, = 0, H, = 0, E, H, = 0, 

Ei Ej = - -E'^ gij, 
1 



Hi Hj - - -H gij. 



(6) 
(7) 

(8) 



The energy-momentum tensor associated with 
Maxwell Lagrangian is given by 



Ff_iu — F^ a F ^ ^ ^ 3^1^ 1 



(9) 



in which F = F^^ F''" ^ 2{H'^ - E'^). Using the above 
average values it follows that Eq. reduces to a perfect 
fluid conflguration with energy density p-y and pressure 
as 



where 



Tf_,v = {p-t + P-t) Vu - p~f g^u, 



Pi = - 2 + 



(10) 



(11) 



The fact that both the energy density and the pressure 
are positive definite for all time yields, using the Ray- 
chaudhuri Eq. (|^), the singular nature of FRW universes. 
Thus Einstein equations for the above energy-momentum 
configuration yield llj] 



where Ao is an arbitrary constant. 



(12) 



Tfiij — —4 Lp Ffj_^Fav + {GLc — L) g^^ 



(14) 



in which Lp represents the partial derivative of the La- 
grangian with respect to the invariant F and simil arly 
for the invariant G. In the linear case, expression (nj) 
reduces to the usual form (||). 

Since we are interested mainly in the analysis of the 
behavior of this system in the early universe, where mat- 
ter should be identified with a primordial plasma Jl5| , p^ , 
we are led to limit our considerations to the case in 
which only the average of the squared magnetic field H^ 
survives |l^,|l^,|l8| . This is formally equivalent to put 
i?^ = in Eq. (M), and physically means to neglect bulk 
viscosity terms in the electric conductivity of the primor- 
dial plasma. 

The homogeneous Lagrangian ( p^ requires some spa- 
tial averages over large scales, as given by Eqs. (§|)"(||). If 
one intends to make similar calculations on smaller scales 
then either more involved non homogeneous Lagrangians 
should be used or some additional magnetohydrodynam- 
ical effec t ||l9| ] should be devised in order to achieve cor- 
relation [pO[ at the desired scale. Since the average pro- 
cedure is independent of the equations of the electro- 
magnetic field we can use the above formulas (^)-(||) to 
arrive at a counterpart of expression ( p^ for the non- 
Maxwellian case. The average energy-momentum tensor 
is identified as a perfect fluid (|lO|) with modified expres- 
sions for the energy density and pressure p^ as 



p^^^H^{l~8aH^ 



Pj^-H^{l- 40aH^). 



Inserting expressions (h5|)-(16) in Eq. (H) yields 



H-^° 



(15) 
(16) 

(17) 



where Ho is a constant. With this resuh, a similar pro- 
cedure appUed to Eq. (||) leads to 



A' 



kHl 



The solution of Eq. (18) is implicitly given as 



(18) 



ci = ± 



"6?" 



(19) 



As far as the right-hand side of Eq. (r8|) must not be 

negative it follows that, regardless of the value of e, for where A{0) = Aq. The linear case ( |l^ ) can be achieved 



a > the scale factor A{t) cannot be arbitrarily small. 
A closed form of Eq. (n3) for e ±1 can be derived as 



from Eq. ( |19D by setting a = 0. 



z — Xi 



(xi - xs)£ (^arcsin ^ y _ 



x-iT ( arcsm 



- Xi 

where xi, X2, are the three roots of the equation SakH^ — kH^x + 3ex^ = 0, and 



--A^it) 



(20) 



T{x, k) 



v/(1-z2)(1-k2^2) 



dz, £{x, k) = 



VT^z 



dz 



are the elliptic functions of the first and of second kind, respectively (see expressions 8.111.2 and 8.111.3 in Ref. pl[). 
The behavior of A{t) for e = ±1 is displayed in the Fig. |l|. 



A(t) 
0.1 
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FIG. 1. Plots from Eq. (|T|). We set A{1) = 1, kH^ = 12, and aH^ = (0; ±1.25 x lO"*) as illustrative values. 



For the Euclidean section, by suitably choosing the ori- 
gin of time, expression ( p^ ) can be solved as 



Expression (21) is singular for a < 0, as there exist a 



= Ho\l^ikc^t^ + Ua). 



(21) 



time t = \J —Ylajkf?- for which A(t) is arbitrarily small. 
Otherwise, for a > we recognize that at i = the radius 
of the universe attains a minimum value A„ 
given from 



which is 



^From Eq. (|T^, the average strength of the magnetic 
field H evolves in time as 



42 



HoVWa. 



(23) 



H' 



1 



2 fcc2t2 + 12a' 



(22) 



Therefore, the actual value of Amm depends on Ho, which 
turns out to be the unique free parameter of the present 



model. The energy density p-y given by Eq. (15) reaches 
its maximum value Pmax = l/64a at the instant t = tc, 
where 



tr 



1 I2a 



(24) 



c V k 

For smaller values of t the energy density decreases, van- 



ishing at t = 0, while the pr essure becomes negative. 
Only for times t ^lQ^JaJk(? the nonlinear effects are 
relevant for cosmological solution of the normalized scale 
factor, as shown in Fig. ^. Indeed, solution ( ^ ) fits the 
standard expression (12) of the Maxwell case at the limit 
of large times. 
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FIG. 2. On the left panel: time dependence of the electromagnetic energy density p., and pressure p-y. pmax = l/64a and 
is given by Eq. (^). On the right panel: nonsingular behavior of the scale factor A{t). Amin and tc are given from Eqs. 
!3|),(|2^. The corresponding classical expression is shown (dashed line) for comparison, with Ao — Amin- 



The energy-momentum tensor ( p^ is not trace-free for 
a 7^ 0. Thus, the equation of state Pj = p^{p^) is no 
longer given by the Maxwellian value; it has instead a 
quintessential-like term ]2^ which is proportional to the 
constant a. That is 



■P7 



16 



(25) 



Equation ( p5|) can also be written in the form 

P7 = ^ Pt - ^ {(1 - 32apT,) 

+ [1 - 2Q{t - tc)Wl~Map-,] , 



(26) 



where Q{z) it the Hcaviside step function. The right- 
hand side of Eq. ( p6|) behaves as (1 — 64ap^)/3^/3 for 
t > tc the Maxwell limit ap^ <^ 1. 

The maximum of the temperature corresponding to 
t = tc is given by 



T„ 



1/4 



,24q;(T/ 

where a is the Stefan-Boltzmann constant. 



(27) 



IV. CONCLUSIONS 



The consequences of the minimal coupling of gravity 
with second order nonlinear electrodynamics (|l3| ) were 
examined. /^From the cosmological point of view, the 
proposed modification is relevant only in the primeval 
era of the universe. Indeed, the class of theories a > 
leads to nonsingular solutions for which the scale factor 
A(t) attains a minimum value. The regularity of the cos- 
mological solution ( ^ll) is to be attributed to the fact 
that, for t < y^24:a/kc^, the quantity p + 3p becomes 
negative. 
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APPENDIX A: ULTRARELATIVISTIC MATTER 
CONTRIBUTION 

Beside photons there are plenty of other particles, and 
physics of the early universe deals with various sort of 
matter. In the standard framework they are treated in 
terms of a fluid with energy density p^, which satisfies 
an ultrarelativistic equation of state Pi, = pii /3. Adding 
the contribution of this kind of matter to the average 
energy-momentum tensor of the photons it follows that 
Pi, = KA~^, where K is an arbitrary positive constant. 
This result allows us to treat such extra matter as noth- 
ing but a reparametrization of the constants Ho and a 
as = + 2K and a = aH^/{H^ + 2Kf . The net 
effect of this is just to diminish the value of Amin as 

1/4 

(Al) 



+ 2K 

Therefore, it turns out that the phenomenon of reversing 
the sign of the expansion factor 3 A/ A due to the high neg- 
ative pressure of the photons is not essentially modified 
by the ultrarelativistic gas. Only an exotic fluid possess- 
ing energy density PexoUc ^ ^" with n < — 8 could be 
able to modify the above result. However, this situation 
seems to be a very unrealistic case. 



[1] E. W. Kolb and M. S. Turner, The Early Universe 
(Addison-Wesley, Redwood City, CA, 1990). 

[2] R. Brandenberger, in Proceedings of the VIII Brazilian 
School of Cosmology and Gravitation, edited by M. Nov- 
ello (Editions Frontieres, Singapore, 1996). 

[3] L. Kofman, A. Linde, and A. A. Starobinsky Phys. Rev. 
D 56, 3258 (1997). 

[4] W. de Sitter, Proc. K. Ned. Akad. Wet. 19, 1217 (1917). 

[5] M. Novello and J. M. Salim, Phys. Rev. D 20, 377 (1979); 
A. Saa, E. Gunzig, L. Br enig, V. Faraoni , T. M. Rocha 
Filho, and A. Figueiredo, |gr-qc/0012105| (2000). 

[6] V. Mukhanov and R. Brandenberger, Phys. Rev. Lett. 68, 
1969 (1992). See also R. Brandenberger, V. Mukhanov, 



[7: 

[9 
[10 

[11 

[12 
[13 
[14 

[15 

[16 
[17 
[18 
[19 

[2o; 

[21 

[22; 



and A. Sornborger, Phys. Rev. D 48, 1629 (1993); R. 
Moessner and M. Trodden, ibid 51, 2801 (1995). 
M. Novello, L. A. R. Oliveira, J. M. Salim and E. Elbaz, 
Int. J. Mod. Phys. A 1, 641 (1993). 

G. L. Murphy, Phys. Rev. D 8, 4231 (1973); J. M. Salim 
and H. R de Olivera, Acta Phys. Pol. B 19, 649 (1988). 
R. Garcia-Salcedo and N. Breton, Int. J. Mod. Phys. A 
15 (27), 4341 (2000). 

R. Klippert, V. A. De Lorenci, M. Novello, and J. M. 
Salim, Phys. Le tt. B 472, 27 (2000); G. Veneziano, 
hep-th/00(D2094| 2000. 



S. W. Hawking and G. F. R. Ellis, The Large Scale Struc- 
ture of Spacetime (Cambridge University Press, Cam- 
bridge, England, 1973); R. M. Wald, General Relativity 
(Univ. Chicago Press, Chicago, 1984). 
M. Novello, in Proceedings of the II Brazilian School of 
Cosmology and Gravitation, edited by M. Novello (J. San- 
som & Cia., Rio de Janeiro, 1980) (in Portuguese). 
R. C. Tolman and P. Ehrenfest, Phys. Rev. 36, 1791 
(1930); M. Hindmarsh and A. Everett, Phys. Rev. D 58, 
103505 (1998). 

H. P. Robertson, Rev. Mod. Phys. 5, 62 (1933); D. Ed- 
wards, Astrophys. Space Sci. 24, 563 (1973); R. Coquer- 
aux and A. Grossmann, Ann. Phys. (N.Y.) 143, 296 
(1982); M. Dabrowski and J. Stelmach, Astron. J. 92, 
1272 (1986). 

T. Tajima, S. Cable, K. Shibata, and R. M. Kulsrud, 

Astrophys. J. 390, 309 (1992); M. Giovannini and M. 

Shaposhnikov, Phys. Rev. D 57, 2186 (1998). 

A. Campos and B. L. Hu, Phys. Rev. D 58, 125021 

(1998). 

G. Dunne and T. Hall, Phys. Rev. D 58, 105022 (1998); 
G. Dunne, Int. J. Mod. Phys. A 12, 1143 (1997). 
M. Joyce and M. Shaposhnikov, Phys. Rev. Lett. 79, 
1193 (1997). 

C. Thompson and O. Blaes, Phys. Rev. D 57, 3219 
(1998); K. Subramanian and J. D. Barrow, ihid 58, 
883502 (1998). 

K. Jedamzik, V. Jatalinic, and A. V. Olinto, Phys. Rev. 
D 57, 3264 (1998). 

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, 
Series, and Products, (Academic, London, 1965). 

R. R. Caldwell, R. Dare, and P. J. Steinhardt, Phys. Rev. 
Lett. 80, 1582 (1998). 



